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ABSTRACT 

We present an analysis of archival X-ray observations of the Type IIL supernova SN 1979C. 
We find that its X-ray luminosity is remarkably constant at (6.5 ± 0.1) x 10 38 erg s _1 over a 
period of 12 years between 1995 and 2007. The high and steady luminosity is considered as 
possible evidence for a stellar-mass (~ 5-10M Q ) black hole accreting material from either a 
supernova fallback disk or from a binary companion, or possibly from emission from a central 
pulsar wind nebula. We find that the bright and steady X-ray light curve is not consistent with 
either a model for a supernova powered by magnetic braking of a rapidly rotating magnetar, or 
a model where the blast wave is expanding into a dense circumstellar wind. 

Subject headings: supernovae: individual(SN 1979C) -supernova remnants 



1. Introduction 



Kasen & Bildstenl (120091) and IWoosleyl (2009) have recently suggested that a class of Type IIL super- 



novae (SNe) may be powered by the birth of a magnet ar during the SN event. Th ese SNe include SN 1961F, 
SN 1979C, SN 1980K, and SN 1985L. In particular, kasen & Bildstenl J2009h argue that for a supernova 
with ejecta mass M e j = 5M that forms a magnetar with an initial period of Pi = 10 ms, the supernova 
can reach peak luminosities consistent with what was observed in these events. They also show that even 
brighter (Mb = -21) events can occur, such as the ultra-bright Type IIL SNe 2005ap and 2008es. 



Motivated by [Kasen & Bildstenl (120091) and IWoosleyl (120091) . we compared the time evolution of the 
X-ray luminosity of SN 1979C to the form expected from a magnetar-powered SN. In § 2, we discuss 
SN 1979C and examine its X-ray luminosity. We find that contrary to the expectation of these papers, 
the evolution of this SN is not consistent with the magnetar model. Additionally, we show that the X-ray 
emission spectrum exhibits evidence for hard X-ray emission, possibly originating from a central source. 
We show that the X-ray luminosity has been steady over the lifetime of the SN evolution and argue that this 
is due to an accreting stellar-mass black hole remnant at the center of SN 1979C. 



2. SN 1979C 



The Type IIL SN 1979C, discovered on April 19 197 9 by G. Johnson dMattei et al.lll979l) and its host 



galaxy NGC 4321 ( M100, at a distanc e of 15.2 ± 0.1 Mpc ; 



served in the radio (We iler et all 19861 : 



Montes et al. 



2000; 



Battel & Bietenholz 



Freedma n et all200ll) have been extensively ob- 



2003, 2008; Marcaide et al. 
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2001 



Immleretal 



2009) . optical jjFesen & Matonick!ll993l : IMilisavljevic et al.l l2009) and X-ray dlmmler et all 1 19981 : iKaaret 



2005) bands. At optical wavelengths, the Ha flux has decreased by ~ 35% between 
19 93 and 2008, while forbid den line emission from lines such as [O III] has increased ~ 50% (see Figure 3 



of IMilisavljevic et al.ll2009D . The increase in the forbidden line emission is evidence of shock heating of 
the ejecta by the reverse shock. Similar to the Ha emission, the r adio emission also shows a decline of a 
factor of ~ 5 over the lifetime of the SN (|B artel & Bietenh olz 2008), though early observations provided ev- 
idence for quasi-periodic oscillations that may be the result of a modulation in th e progenitor's circumstellar 
environment by a binary companion (We iler et all 19921 : ISchwarz & Pringldl 1 9961) . 



Recently, IB artel & Bietenholz! (120081 ) noted that the in conjunction with a decrease in the radio lumi- 
nosity the spectrum has begun to flatten. They note that the flattening spectrum would be expected for a 
supernova with emission from a compact remnant that is beginning to appear in the shell as the expanding 
remnant becomes increasingly transparent. They point out, however, that the flattened spectrum could also 
be due to synchrotron emission from the shell as electrons are accelerated to GeV energies. We note that 
in this scenario, it would be reas onable to expec t an i ncreasing hard X-ray flux as electrons are accelerated 
to TeV energies. More recently, iMarcaide et al.l (|2009r) presented a revised analysis of the radio expansion 
of SN 1979C and found that the expanding blastwave is in near free expansion with m = 0.91±0.09 for 
R (x t m . This seems contrary to X-ray observat ions which suggest that the blastwave has been strongly 
decelerated by a dense circumstellar interaction dlmmler et al.1 120051) . While the optical and radio emis- 
sion have showed evolution over ~ 25 years, X-ray observations of SN 1979C show a markedly different 
behavior. As seen in Table[T] the X-ray luminosity has remained remarkably constant over its 30 yr lifetime. 



Kasen & Bildstenl (|2009l ) suggest that Type IIL SNe such as SN 1979C may have been magnetar 



powered. In the context of a magnetic dipole model, one can write the spindown luminosity as L p = 
5 x lO^Bf^P^ 4 erg s _1 for a magnetar with magnetic field B14 (in units of 10 14 Gauss) and period Pi (in 
units of 10 ms), or as 



L p = 2 x 10 42 5 2 4 (t/yr)- 2 erg s" 1 . (1) 

The observed X-ray luminosity, Lx, can be expressed as some fraction, f x , of the magnetar spindown 
luminosity, i.e. Lx = f x L p . In FigureQ]we plot the evolution of the X-ray luminosity for a dipole magnetar- 
powered SN assuming that 1, 3, or 10% of the spindown power is transferred into X-rays, with E>n = 1 and 
P i = 1 . We also indicate the X-ray luminosity (from Table [T]) as observed by the ROSAT, XMM-Newton, 
Chandra, and Swift observatories, as well as the early upper limit from the Einstein observatory HRI, where 
the model predicted X-ray luminosity is more than two orders of magnitude above the Einstein upper limit. 
As seen in Figure Q] and also in Table [T] the observed X-ray luminositjQ is remarkably constant over time, 
withLx = (6.5±0.1)xl0 38 erg s -1 . 

Most importantly, even if the magnetar powered the SN only at early times, it would have likely resulted 



'to convert from the observed count rates to source fluxes and luminosities, we assume an absorbed thermal bremsstrahlung 
model with n H = 2.4 x 10 20 cirT 2 and kT = 0.5 keV. 
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in a clear detection by the Einstein observatory. The above equation applies to a magnetic dipole spindown. 
More generally, the magnetar luminosity can be written as 

(Z-l) 

Lp a (i+t/t P r (2) 

where t p is the initial spindown time and I = 2 for magnetic dipole spindown. At late times, L p oc t , 
and for all physical values of I the luminosity decreases with time. If the X-ray luminosity is some fraction 
of the spindown luminosity, the fraction will need to increase with time, potentially up to unphysical values 
greater than unity, in order for the observed X-ray luminosity to remain constant. 



2.1. Spectral Modeling 



The combined Chandra ACIS-S spectrum of SN 1979C is shown in Figure |2l llmmler et al.1 (12005!) fit 
the XMM-Newton observation to a two component thermal plasma model (MEKAL model in XSPEC0), with 
kTim,, = 0.77!q 19 keV and kThigh = 2.31^ 'gg keV. They did not find any evidence for strong absorp tion, 



- high 



with a fitted hydrogen column density Nh ~ 2.5 x 10 20 cm -2 (see also iDickey & Lockmanlll990l) . and 
attribute the low temperature component to recently shocked ejecta, and the high temperature component to 
emission from swept-up circumstellar material (CSM). In their model, the shocked circumstellar emission 
arises from a dense wind, with den sities of < 10 5 cm~ 3 out to radii of 0.13 pc, yielding column densities 



< 10 22 cm" 2 . 



Immler et all (|2005h argue that VLBI imaging at t = 22 yr show no evidence for clumping in 



the CSM, indicating that the CSM is remarkably uniform. If this were the case, then both the XMM-Newton 
and Chandra ACIS-S spectrum would show strong absorption at low energies. We fitted the Chandra ACIS- 
S spectrum to a two component thermal plasma model, and found that the spectrum can be described by two 



thermal plasmas with temperatures kT = 0.751q 2 and 1.86^q'°2 keV, which shows no evidence for the 
level of absorption expected fro m the above propertie s of a hypothetical dense smooth CSM. We note that 
our fitted parameters and those of llmmler et all (|2005r) are consistent. 



0.88 



Recently, IBietenholz et al.1 (|2010T ) reported on VLBI detection of a central source in SN 1986J. They 
argue that the small proper motion of the source (500 ± 1500 km s -1 ) provides evidence that the source 
is either a central neutron star or black hole, but they do point out that it might also be a dense CSM 
condensate seen along the line of sight, as similar clum ped regions app ear in the northeastern region of the 
SN 1986J shell. Additionally, as previo usly pointed ou t, Marcaide et al.1 (120091) showed that the blastwave of 



SN 1979C is still in free expansion, and lBartel & Bietenhola ([2008) argued that the flattening radio spectrum 
could be evidence for emission from a central compact remnant. We thus chose to fit the X-ray spectrum 
of SN 1979C as a combination of shocked plasma, either arising from shocked CSM or ejecta, along with a 
component associated with a central source. The resultant fits are shown in Figure|2l In both cases, we found 
that the thermal component is well described by a plasma with kT = 1.1 ^q'^ keV in collisional ionization 



2 http://heasarc.nasa.gov/xanadu/xspec/ 
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equilibrium. We first modeled the hard emission as a power-law (bottom panel of Fig. [2]), and found a 
power-law index for the flux per unit energy T = 2.2+q"|. We also fit the data to a model for an accreting 
central black hole, using the kerrbb model in XSPEC, shown in the top panel of Figure [2 We found that 
the hard component of the spectrum can be well fit for a black hole mass of ~ (5.2 ± 0.8)M Q . We do note, 
however, that while the X-ray spectrum is well described by a spectrum from an accreting black hole, we 
do concede that a powerlaw fit such as that shown in the bottom panel of Figure |2] is still consistent with 
emission from either a central compact source such as a pulsar wind nebula, or from synchrotron emission 
arising from shock accelerated electrons. 



2.2. A Black Hole Remnant? 

While a two component thermal model can describe the X-ray spectrum for SN 1979C, with the hard 
component having a kT > 2 keV, as seen in Figure [2 the data are equally well described by a thermal 
component, either from shocked CSM or shocked ejecta, and a hard spectral component, possibly from 
a central source. Ilmmler et al.l (120051 ) modeled the X-ray emission as arising from the expansion of the 
blastwave into a dense circumstellar wind, with values ranging from 10 4 - 10 7 cm -3 out to radii of 4x 10 17 
cm. Ahead of th e blastwave, the circ umstellar density in their model is still 10 3 cm -3 at a distance of 1 pc 
(c.f. Figure 6 of Ilmmler et al.ll2005r) . which would add an absorbing column density of ~ 3 
well in excess of the measured absorption in the spectrum of SN 1979C. 



10 21 cm" 2 , 



Additionally, in models for supernova remnant (SNR) evolution where the blast wave is expanding into 
a wind, the X-ray (and optical and radio) emissivity from the expanding blast wave should decrease with in- 
creasing blast wave radius. This is because the blast wave encounters the densest material at small radii, and 
sweeps over less dense material as it expands to a larger volume. At the late times of interest here, the ejecta 



is optically thin (IChevalier & Franssonlll994l) . The free-free emission measure scales quadratically with the 
particle number density n and linearly with the volume V oc r 3 , so that for the typical n oc r~ 2 radial profile 
of a wind, Lx oc n 2 V oc r _1 decreases with time. Other SN remnants (SN Rs) do show evidence for a strong 
intera ction with a circumstellar medium (c.f., SN 19941 and SN 1993J: Ilmmler et al.ll2002k IChandra et al. 
l2009h . SN 1993J was observed with Chandra in 2000 and again in 2008. Analysis of these data show that 
the X-ray luminosity has dropped from ~ 6 x 10 38 erg s _1 in 2000 to ~ 1.5 x 10 38 erg s _1 in 2008. This 
sharp drop in X-ray luminosity is inconsistent with the near steady X-ray emission from SN 1979C. 



A plausible explanation for the constant luminosity would be to associate it with Eddington-limited 
accretion onto a central compact object of mass M x . Based on the Eddington luminosity value, 



L Edd = 1.4 x 10 



.88 



M 



x \ 1 

erg s , 



(3) 



with a modest bolomet r ic correction (e.g. Lfj n = 1 .6xl0 37 , Lnv = 9xl0 36 and Li.qghz = 1.6xl0 36 



-i. 



Immler et al.1 120051 : iBartel & Bietenhola 120031) . this implies that the accreting object has a mass 



erg s 

of 5-lOM , which is within the typical range associated with stellar-mass black holes (|McClintock et al 
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2009). The hard component of the X-ray spectrum of the source shown in Figure [2] is consiste nt with the 



spect rum of an Eddington-limited black hole in an X-ray binary such as LMC X-3 (see Fig. 1 in iDavis et al. 



2006), and our spectral fits to the data imply a black hole mass of > 5M^. It is p lausible to expect that 
a black hole forms in a type ILL SNe. According to Figure 2 of LHeger et al.1 (120031) . so lar metallicity stars 
with progenitor masses ~ 25M Q can follow a track that leads to black hole formation. iHeger et all (120031) 
also point out that 25% of low mass (M ~ 20M©) progenitors with solar me tallicity will p roduce black 
hole remnants by fallback. Since ~ 3% of core-collapse SNe are Type IIL (jSmarttl 12009). this implies 
that some Type IIL's will form a black hole remnant, consistent w ith predictions f or theoretica l black hole 
mass distributions from core-collapse SNe (IFryer & Kalogerall200ll) . Interestingly, Heg er et all ([2003) also 
suggest that Ty pe IIL/b SNe are produced in binarie s, which would avoid the need for a fallback accretion 
disk in the SN (IPerna et alj|2000l : IWang et al.ll2006T) as the binary companion would provide the accreting 
material. There has been some evidence for a binary companion to the progenitor of SN 1979C, seen as a 
modulation of the early time radio light curve. 

The existence of an accreting black hole remnant at the center of SN 1979C is consistent with the 
upper limit on the early X-ray luminosity established by the Einstein Observatory. For supernova models 
with 2-5M Q of ejecta and explosion energies of 2 x 10 51 ergs, the expanding ejecta bubble becomes optically 
thin to keV X-rays on a timescale of 15 yrs. The upper limit on the Einstein observation may reflect that 
early phase of the SNR evolution when the ejecta is optically thick or the time delay required to establish a 
stable accretion fl ow around the black hole remnant, which i s calculated to be about twice the recombination 
timescale, ~ 1 yr ( Zampieri et al. 1998 : Balberg et allbood) . 



Milisavljevic et al.1 (|2009l) note the presence of a Wolf-Rayet "bump" in the optical emission spectrum 
of SN 1979C which may be associated with the progenitor . They mention tha t several metal-rich H II 
regions in Ml 00 have been observed to possess WR stars (Ivan Dyk et al.lll999i : IPindao et al.ll2002l) . and 
HST imaging shows several young blue stars of age 4-6 Myr in the vicinity of SN 1979C. Photometry of the 
clusters' stellar population yields an estimate of the progenitor's mass to be (18)±3 M . This is consistent 
with the expected progenitors for Type IIL SN and with the required progenitor mass to form a black hole 
during the SN explosion, possibly from a fallback disk. 

We have also analyzed Chandra and ROSAT observations of other two other Type IIL SNe, SN 1980K 
and SN 1985L. We have found that the X-ray luminosity of SN 1979C is more than an order of magnitude 
larger than that of SN 1980K (L x « 5.5 xlO 37 erg s" 1 ) and SN 1985L (L x < 4x 10 37 erg s" 1 ). Inter- 
estingly, the latter luminosities from these Type IIL S Ne are more consistent with a pulsar or igin in young 
SNe, such as SN 1968D, SN 1941C, and SN 1959D JPerna et al.lbooil : Isoria & Pernall2008h . in line with 
Heger et all (120031) who point out that not all Type IIL's will form a black hole remnant. 



3. Conclusions 



Our analysis of archival X-ray observations of SN 1979C indicate that the X-ray luminosity has been 
remarkably steady. We find that the X-ray light curve is not consistent with either a model for a supernova 
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powered by a magnetar or a model where the X-ray emission arises from a blast wave expanding into a 
dense circumstellar wind. In the latter case, the observed decline in the optical and radio bands should be 
accompanied by a decline in the X-ray flux. The X-ray spectrum for SN 1979C can be modeled by a com- 
bination of a thermal X-ray component and emission from an accreting black hole with a mass ~ 5M . The 
accreting material likely originates from either a fallback disk after the supernova, or possibly f rom material 



accre t ed from a binary compan ion which is suggested by the radio light curve at early times (|Weiler et al 



19921 : ISchwarz & Pring le 1996). Finally, we note that the format ion of a black hole in SN 1979 C might 



have also imprinted a feature in the optical supernova light curve. lYoung. Smith. &Johnsonl (120051) recently 
showed that the supernova light curve can be fit by a two component model which includes a GRB afterglow 
followed by supernova ejecta. They argued that the GRB optical afterglow is produced when a jet, from the 
formation of the central ~ 2M Q black hole, penetrates through the stellar envelope. 

We note that SN 1979C appears to be the first example of a historic supernova where there is possible 
evidence for a black hole remnant. A survey of Type IIL SNe observed at late times could reveal the 
existence of other accreting black hole remnants and constrain the statistics of black hole formation in core- 
collapse SNe. We note that a deep Chandra or XMM-Newton observation of SN 1979C could settle the 
issue, as it would allow for both a detailed spectral analysis of the emitted spectrum as well as a search for 
short term variations in the X-ray light curve, which could be seen as evidence for ongoing accretion. 
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Table 1. X-ray observations of SN 1979C 



At 


Count Rate 






Mission 


yr 


10" 4 cps 


10~ erg cm -2 s^ 1 


10 38 erg s^ 1 




0.7 


< 3.0 


<2.3 


< 6.3 


Einstein (HRI) 


16.2 


6.7±0.7 


3.0±0.3 


8.2±0.9 


ROSAT (HRI) 


20.6 


42.±2.0 


2.5±0.2 


6.9±0.6 


Chandra (ACIS-S) 


22.7 C 




2.3±0.3 


6.3±0.7 


XMM-Newton (MOS) 


26.5 d 


8.0±0.9 


2.3±0.3 


6.3±0.7 


Swift (XRT) 


26.9 


40.±0.8 


2.4±0.2 


6.6±0.5 


Chandra (ACIS-S) 


28.0 


43. ±0.3 


2.6±0.2 


7.0±0.5 


Chandra (ACIS-S) 



a Unabsorbed 0.3-2.0 keV flux assuming a 0.5 keV thermal Bremsstrahlung model 
with Galactic hydrogen column density of 2.4 x 10 20 cm -2 . 



b uminosity assumes a distance to M100 of 15.2±0.1 Mpc (IFreedman et al.ll2001l) . 



c Flux value taken from llmmler et al.l (120051) . 

d Values are from Swift XRT observations taken between 2005 and 2006, with the 
first observation performed on 2005-11-13. 
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Fig. 1. — Evolution of the X-ray luminosity of SN 1979C. The dotted, dashed, and long-dashed diagonal 
lines correspond to the evolution of the SN luminosity if it were powered by the magnetic dipole spindown 
of a central magnetar for various X-ray efficiencies, f x . The data points correspond to the 0.3-2.0 keV X-ray 
luminosities with 1-cr uncertainties, assuming a distance of 15.2 Mpc. The upper limit from the Einstein 
observatory is depicted as an arrow. The heavy solid line corresponds to our best fit value of a nearly constant 
X-ray luminosity. 
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Fig. 2. — The unfolded Chandra ACIS-S spectrum of SN 1979C. Panel (a) shows the data with its l-a error 
bars (crosses), which is fitted as a combination of a thermal plasma component (blue long-dashed line), and 
emission from a relativistic accretion disk around a Kerr black hole (red short-dashed line). For the thermal 
component, we infer a temperature of kT = 1.2 keV, and for the disk component, we infer a black hole 
mass of Mbh = 5.2M Q . Panel (b) shows the same thermal plasma model as in panel (a) together with a 
power-law fit (red short-dashed line) with a spectral slope for the flux per unit energy Fe oc E~ r of T = 2. 1. 
In both panels, the fitted absorption is consistent with Galactic foreground of a hydrogen column density, 
iV H = 2.5 x 10 20 cm" 2 . 



